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Abstract 

We present results from Monte Carlo simulations of a three dimensional fermion field theory 
which can be derived from a model of graphene in which electrons interact via a screened 
Coulomb potential. For our simulations we employ lattice gauge theory methods used in ele- 
mentary particle physics. We show that the theory undergoes a second order phase transition 
and we provide estimates for the critical exponents. The estimated value of the physical critical 
coupling implies that graphene in vacuum is an insulator. We also present the first results for the 
quasiparticle dispersion relation. 

Key words: graphite, electrical conductivity, simulation 



1. Introduction 

There has been considerable recent interest in graphene (single layer graphite) sparked by its 
discovery and subsequent experimental study [1]. The remarkable properties of graphene have 
been suspected for many years [2]. In brief, for a carbon monolayer having one mobile electron 
per atom, a simple tight-binding model shows that the spectrum of low-energy excitations ex- 
hibits a linear dispersion relation centred on zeros located at the six corners of the first Brillouin 
zone (see e.g. (3[]). Using a linear transformation among the fields at two independent zeros it is 
possible to recast the Hamiltonian in Dirac form with Nf = 2 flavors of four-component spinor 
d) - the counting of the massless degrees of freedom coming from two carbon atoms per unit cell 
times two zeros per zone times two physical spin components per electron. Electron propagation 
in the graphene layer is thus relativistic, albeit at a Fermi velocity Vp ~ 10 6 m/s. The implications 
for the high mobility of the resulting charge carriers is the source of the current excitement. The 
electronic properties of graphene are discussed in detail in a recent review paper |4j]. 

Although the semimetallic properties of graphene on a substrate are well-known, it remains 
unclear whether graphene in vacuum undergoes a phase transition from the semimetal phase to an 
insulating phase. A sufficiently strong Coulomb interaction may lead to a condensation of elec- 
tron and hole pairs thus turning the conducting phase into a gapped insulator at a critical value of 
fermion flavors N/ c . Recent estimates of Nf c have been obtained by: self-consistent solutions of 
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Schwinger-Dyson equations yielding N/ c = 2.55 JH; renormalization group treatment of radia- 
tively induced four-fermion interactions yielding N/ c — 2.03 |@1; and lattice simulations yielding 
Nf c = 4.8(2) |7]. Furtermore, a lattice simulation of graphene based on a four dimensional gauge 
invariant action with fermions propagating in two spatial dimensions predicted that freely sus- 
pended graphene is an insulator 18j,|9|]. The results of lattice simulations of a Thirring-like model 
pertinent to graphene presented in this paper are consistent with this observation. 



2. The Model 



While the above considerations apply quite generally, a realistic model of graphene must 
incorporate interactions between charge carriers. One such model, due to Son II 1011 . has Nj mass- 
less fermion flavors propagating in the plane, but interacting via an instantaneous 3d Coulomb 
interaction. In Euclidean metric and static gauge 8qAq = the action reads 
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where e is the electron charge, V = Aq is the electrostatic potential, and the 4x4 Dirac matrices 
satisfy {y^, y v ) = 2<5 A , V , yU = 0, . . . , 3. In our notation x is a vector in the 2d plane while the index i 
runs over all three spatial directions. In the large-A^/ limit the resummed V propagator becomes 
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where p 2 = Pq + v^p] 2 - Dq corresponds to the classical propagation of V and n comes from 
a vacuum polarisation fermion-antifermion loop. In either the strong coupling or large-A^/ lim- 
its D\(p) is thus dominated by quantum corrections - the relative importance of the original 
Coulomb interaction being governed by a parameter A = |n/Do| Po =o- Chiral symmetry breaking 
due to the spontaneous condensation of particle - hole pairs is signalled by an order parameter 
(i^i/r) + 0. Physically the most important outcome is the generation of a gap in the fermion 
spectrum, implying the model is an insulator. 

The physics of eq. (JTJ is very reminiscent of the three dimensional Thirring model, which is 
analytically tractable at large-A^/, but exhibits spontaneous chiral symmetry breaking at small Nf 
and strong coupling [11, 11211 . Arguably the Thirring model is the simplest field theory of fermions 
requiring a computational solution. The location of the phase transition at Nf = Nf c in the strong 
coupling limit has recently been determined by lattice simulations to be Nf c = 6.6(1) [13]. The 
apparent similarity of the two systems has led us to propose a Thirring-like model pertinent to 
graphene H, with Lagrangian 
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As for eq. ([T} we assume a large-A^/ limit to estimate the dominant vacuum polarisation correc- 
tion; the resultant propagator for V is 
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Figure 1: (if/ijj) versus g 2 with fits to eq. j5J- 



In the strong-coupling or large-A^/ limits, D2 coincides with D\, implying that the fermion inter- 
actions are equivalent. We thus expect predictions made with the model eq. (0), and in particular 
critical behavior such as the value of Nf c = 4.8(2) [7], to be generally valid for Son's model 
eq. ([T} in the limit of large A. 



3. Lattice Simulation Results 

In this section we present results from simulations with Nf = 2 which corresponds to real 
graphene. Details for the lattice action and the algorithm employed for our simulations can be 
found in 0|. The anisotropic nature of the model dynamics results in the systematic effects 
due to finite temporal extend L, being more important than those due to finite spatial extend L s . 
By comparing the values of the chiral condensate from simulations on different lattice sizes we 
found that the finite size effects on 24 2 x 48 are within statistical errors for the parameter ranges 
m = 0.005, 0.525 < g' 2 < 0.65 and m = 0.010, 0.025, 0.525 < g~ 2 < 0.70 (where m is the 
fermion bare mass). The chiral condensate as a function of mo and (g~ 2 - g~ 2 ) near a continuous 
phase transition scales according to an equation of state given by 

m = A(g- 2 - g- 1 )®*) 8 -* + B(W) S (5) 

We used eq. © to curve-fit our data and extracted g~ 2 = 0.607(5), 6 = 2.72(7), and p = 0.68(2). 
The data and the fitting curves including the mo = curve are shown in Fig. Q] Our estimates 
of the critical exponents are clearly distinct from their strong coupling values 6 = 5.5(3) and 
P = 0.22(2) J0]. We also note that the exponent 6 is larger than the estimate 5 = 2.26(6) of Drut 
and Lahde 19|]. 

The fermion dynamical mass my is an inverse time-like correlation length of the system. At 
the critical coupling, (jfnlr) and m j obey the following scaling relations for mo ± 0: 

{xj/if/) - aiinj 6 , mf=a2in v . (6) 



3 



0.25 - 



o Chiral condensate 
• Fermion mass 



0.2 - 



0.15 



0.1 - 



*'©-' 



..o-~ 



0.05 - 



0.01 



0.02 
m 



0.03 



0.04 



Figure 2: (if/i//) and tnt fits to eq. (6) at/3 = 0.60. 



By fitting our data at /3 = 0.60 (which is close to the critical coupling g~ 2 = 0.607(5)) to eq. © 
we obtain 6 = 2.81(2) and v, - 0.80(6). This value of 6 is consistent with 6 = 2.72(7) obtained 
from the equation of state. In Fig. [2] we present the data and the fitting functions. 

For systems with anisotropic interactions the correlation lengths in the time and space direc- 
tions may scale with different exponents v, and v s respectively. The dynamic critical exponent 
z = v,/v s describes relative scaling of time and space length scales and governs the quasiparticle 
dispersion relation E ~ p z for p — > 0. In d spacetime dimensions the exponents v, and v s are 
related to /? and 6 via a generalized hyperscaling relation 11411 : 

v, + (d-l)v s =0(6+1) (7) 

From eq. (0 we obtain v s = 0.90(5) and z = 0.89(8). This value of z is close to the strong 
coupling and large-A^/ z ~ 0.8 estimate of Son Jioll . A value of z < 1 implies that the quasiparticle 
excitations are stable, because energy momentum conservation forbids their decay into two or 
more quasiparticles. 

Before we decide whether our model predicts that graphene in vacuum is a conductor or an 
insulator we need an estimate of the physical critical coupling A c . First we calculate the renor- 
malized critical coupling g~%. We must note however, that the vacuum polarisation calculation 
leading to eqs. (fJJ and dU does not go through in quite the same way for the lattice regularised 
model; rather, there is an additive correction which is momentum independent and ultraviolet- 
divergent: 

n la "(p) = n co "i( P ) + g 2 j(m), (8) 



where J(m) comes from incomplete cancellation of a lattice tadpole diagram [11]. This extra 
divergence not present in the continuum treatments can be absorbed by a wavefunction renor- 
malisation of V and a coupling constant renormalisation 

Sk = 1 8 lu v (9) 

1 - g z J(m) 

In the large-Af/ limit we thus expect to find the strong coupling limit of the lattice model at 
g 2 R — > oo implying g 2 — > gi = J~ l (m). For g 2 > gi D ,a "(p) becomes negative, and S i att no 
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Figure 3: Dispersion relation at g 2 = 0.8, 0.9 and mo = 0.005 on a 48 3 lattice with fits to eq. Jilt . 



longer describes a unitary theory. For small-A^/ we identify g 2 im with the value of the coupling 



where the chiral condensate has its maximum value J3l. For Nt = 2 



0.30 A and from 



eq. © we obtain g C R - 3.23(4). Next, we match D}? tt (p) to D\(p) at a reference momentum 
Po = 0, \p\ = f > so that the propagators are equal at a distance of roughly one lattice spacing (for 
shorter distances D 1 "" > D\ and vice-versa) yielding the condition 
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A c = = 2.54(2). 



(10) 



This value of /l e is larger than the value A c « 1.25 calculated for graphene on an Si02 substrate 
fioll . but smaller than the value Aq « 3.4 based on a substrate with unit dielectric constant. 
Our result thus suggests freely-suspended graphene is an insulator. Drut and Lahde obtained 
from simulations of a gauge invariant action of graphene A c — 1 .67(2), which also implies that 
graphene in vacuum is an insulator [8[] . 

A major advantage of our formulation is that it permits a relatively straight-forward analysis 
of the quasiparticle dispersion relation. Here we present our first results from simulations on 
32 2 x 48 and 48 3 lattices. We fitted the fermion energy as a function of momentum p in the x 
direction to the following discretized free field dispersion relation which assumes z = 1 : 



E(p) 



A sinh 1 ( A /sin 2 /;> + M 2 ) 



(ID 



where AM is the quasiparticle dynamical mass m f and A ~ lim p ^o f| is the Fermi velocity. 
Equation (TTTT) fitted the data relatively well, providing no evidence for z i 1 1 . An example of fits 
to data generated with g~ 2 = 0.8, 0.9 and mo = 0.005 on a 48 3 lattice is shown in Fig. [3] In Fig. [4] 
we show that within the accuracy of our analysis the parameter A is independent of both mo and 
g~ 2 with a value roughly A su 0.7 suggesting that the free field value of Vf is renormalized. 
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Figure 4: Values of A versus mo extracted from fits with eq. jilt on 32 2 X 48 and 48 lattices. 



4. Summary and Conclusions 

In this paper we studied a (2+ l)d Thirring-like model pertinent to graphene based on eq. O]). 
We performed Monte Carlo simulations using lattice gauge theory techniques and showed that 
the model undergoes a second order phase transition with critical exponents 6 = 2.72(7), ft = 
0.68(2), and v, = 0.80(6). Using a generalized hyperscaling relation valid for theories with 
anisotropic interactions we obtained v s = 0.90(5) and z = 0.89(8). Our preliminary z < 1 
estimate implies that the quasiparticles are stable based on energy-momentum conservation. The 
value of the physical critical coupling A c - 2.54(2) implies that freely-suspended graphene is an 
insulator in agreement with Jit]. We also presented the first results for the quasiparticle dispersion 
relation and within the accuracy of our simulations we showed that the free field Fermi velocity 
is renormalized. Simulations with higher resolution in momentum on larger spatial lattices are 
currently in progress for a more accurate study of the dispersion relation. 
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